Roughening Transition of Interfaces in Disordered Systems 
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The behavior of interfaces in the presence of both lattice pinning and random field (RF) or random 
bond (RB) disorder is studied using scaling arguments and functional renormalization techniques. 
For the first time we show that there is a continuous disorder driven roughening transition from 
a flat to a rough state for internal interface dimensions 2 < D < 4. The critical exponents are 
calculated in an e-expansion. At the transition the interface shows a superuniversal logarithmic 
roughness for both RF and RB systems. A transition does not exist at the upper critical dimension 
D c — 4. The transition is expected to be observable in systems with dipolar interactions by tuning 
the temperature. 

PACS numbers: 68.35.Ct, 05.20.-y, 68.35.Rh 



The thermal roughening transition (RT) from a flat, 
localized to a rough, delocalized interface in the three di- 
mensional Ising and related lattice models is one of the 
paradigms of condensed matter physics jlj. It can be 
directly observed if a crystal surface undergoes a tran- 
sition from a faceted to a smooth shape At the 
RT temperature Tr the free energy of a step on the sur- 
face vanishes. Apart from its shape also other physi- 
cal properties are influenced by the presence of the RT, 
e.g., the (growth) velocity u of the weakly driven in- 
terface (surface) changes dramatically at the RT from 
u ~ exp(-C//) for T < T R to u ~ / for T > T R . Here / 
denotes the driving force density, which is proportional 
to the magnetic field in the case of Ising magnets and to 
the difference A/i of the chemical potentials of the crystal 
and its melt in the case of crystal growth, respectively. 
The thermal RT was shown to be in the same universality 
class as the metal-insulator transition of a 2D-Coulomb 
gas the Kosterlitz-Thouless transition of a 2D XY- 
model H (both with an inverted temperature axis) and 
the phase transition of the 2D sine-Gordon model ||. In- 
terestingly, Kosterlitz j?]] and Forgacs et al. ||] considered 
the thermal RT for D < 2 interface dimensions and found 
a superuniversal logarithmic roughness at the transition. 
Finally, the thermal RT disappears in D = 1 dimensions: 
1-dimensional interfaces are rough at all finite T M. 

Disorder is an even more efficient source of interface 
roughening and physically important as well |i"o(| . Inter- 
face roughening due to disorder was considered to de- 
termine the lower critical dimension of the random field 
Ising model |TT|JT^ ] and the mobility of domain walls in 
disordered magnets |l3|-|lE|]. It was argued that in the 
presence of disorder interfaces of dimensions D < 2 are 
always rough, but undergo a roughening transition as a 
function of the disorder strength for D > 2 dimensions 
p6[ . Bouchaud and Georges considered explicitly 
the competition between lattice and impurity pinning in 
2 < D < 4 dimensions using a variational calculation. 



Surprisingly they found three phases: a weakly disor- 
dered flat, a glassy rough (GR) and an intermediate flat 
phase (GF) with strong glassy behavior. The transition 
between the GR and the GF phase turned out to be first 
order. For D = 2 only the glassy rough phase was ex- 
pected to survive, but that shadows of the two other 
phases would still be seen in the short scale behavior. 
However, it should be kept in mind that the variational 
calculation is in general an uncontrolled approximation 
which is known to give spurious first order transitions 
p8| and Flory-like exponents as indeed found in . 

It is therefore the aim of the present paper to recon- 
sider this problem by using a functional renormalization 
group calculation in D = 4 — e interface dimensions. It 
turns out that the glassy flat phase found in JItJ is re- 
placed by a crossover region with logarithmic roughness 
around the RT. The model we consider is that of p7| : 



d D 2 



\{Vzf + V L (z(x)) + V fl (z(x),x) 



(1) 



Here z(x) denotes the displacement of the D dimen- 
sional interface from a flat reference configuration, and 
7 is an elastic stiffness constant. For simplicity, the lat- 
tice potential can be considered to be given by Vl(z) = 
—va 2 cos(2irz/ a); higher harmonics could be taken into 
account but turn out to be irrelevant. Below, we will 
measure all transverse lengths in units of a/2ir and there- 
fore put a = 2tt in the following. The random pinning po- 
tential Vr(z(x), x) is assumed to be Gaussian distributed 
with a zero average value Vr{z, x) = and 



V R (z, x)Vr(z', x') = R(z - 2>)&{tl - x'). 



(2) 



The form of R(z) depends on the type of the disorder; in 
particular we have to distinguish between random field 
(RF) and random bond (RB) disorder. For RF's, to be- 
gin with, R{z) ~ — Aq|z| for large z where V Aq denotes 



the strength of the random field |11|. For RB's the bare 
form of R{z) is a smeared out ^-function of width £o and 



1 



strength Ao£q ■ For simplicity we assume that £o is of the 
order a. 

We begin our calculation with some elementary con- 
siderations. For Ao = a step of height 2ir in the oth- 
erwise planar interface of linear dimension L costs an 
energy E stC p — jy/vL ^ 1 . The step itself has a width 
ui s tcp — v^ 1 / 2 . Here and below we omit all prefactors of 
order unity. Switching on RF disorder, we consider the 
energy of a hump of total height z in the interface con- 
sisting of z/2n such steps. This gives for large enough 
lattice pinning v 



E, 



total 



(3) 



The second part represents the interaction energy with 
the RF 0. Minimizing (|) we get z ~ (L/L R )( 2 - D *> 
with L R ~ (v/A ) 1/(2 ~ D) . For D < 2 and L > L R the 
interface is therefore rough even in the presence of the lat- 
tice potential. Since the total step free energy vanishes 
now on scales L ~ L R , the system is described on larger 
scales by an effective elastic Hamiltonian with a stiffness 
constant 7 c ff . A rough estimate for 7e ff follows from bal- 
ancing the elastic and the bare step free energy on the 
scale L Rl which gives 7 c ff ~ jL^/wstcp 16|- On scales 
L ^> L R the hump height scales as z ~ (L/L R ) C ' where 
C = (4 - D)/3 To be meaningful, L R should be 

larger than the bare stepwidth, L R > w s tcp, a condition 



which can be rewritten as v > 



A 



2/(4-D) 



L R can 



then be cast in the more conventional form of a Larkin- 
length L R = (v/v^/^^A- 1 ^ 4 -^. If v -> v c +, 7cff 
changes into 7. For v < v c we can neglect the lattice 
pinning term in (|l]) completely and v has to be replaced 
by v c in the above formulas. 

So far we have considered the case D < 2. For D = 2 
it was shown by Binder |l2|,[l9| that L R becomes expo- 
nentially large, L R ~ exp(Cw/Ao), with C a numerical 
factor, and interfaces are again rough for L > L R . Thus, 
there is no flat phase in D < 2 dimensions. 

For D > 2 we have to expect from the previous consid- 
erations that the interface becomes rough only for v < v Cl 
i.e., that there is a RT as a function of i> (or Ao) at v ~ v c . 
In the rough phase the mean square displacement of the 
interface can be estimated again from the energy 



7~L^ 



z 2 + V^e- z2 / 4 - v/A^i (4 ^ )/2 . (4) 



Here we have included both the elastic and the step 
energy term to describe the crossover from v <C v c to 
v ~ v c . Moreover, we have treated the lattice pinning 
term in a self consistent harmonic approximation to ac- 
count for a situation in which z > 2-7T. In the rough 
phase, the hump height asymptotically scales with the 
roughness exponent Q. However, there is an intermedi- 
ate length scale region in which the second term of (^) 
is larger than the first term and the roughness increases 
only logarithmically, z 2 ~ 2(D — 2)h\(L/ L R ). Finally 
for »»u c the interface is always flat. Similar arguments 
apply to RB systems but the asymptotic roughness can- 
not be obtained from a Flory argument. Below, we will 



show by a detailed renormalization group (RG) calcula- 
tion in D = 4 — e dimensions that the qualitative picture 
obtained so far is correct and that the crossover from log- 
arithmic to power law roughness happens at length scales 
L~£|| ~( Uc - v )-VVi 

As was first shown in [po[ , the proper treatment of 
interface roughening requires a functional RG approach 
due to the occurrence of infinitely many relevant opera- 
tors. A finite lattice pinning potential Vl(z) breaks the 
continuous translational symmetry and thus a functional 
RG treatment of both pinning potentials generates also 
contributions to the elastic stiffness. Balents and Kar- 
dar (2l]] have developed a RG approach in terms of an 
eigenfunction expansion of the interactions for a similar 
case. For the periodic lattice potential, this expansion is 
simply a Fourier expansion. This simplifies the analysis 
of the initial functional form of the RG flow since higher 
harmonics cosmz, m > 1, are strongly irrelevant at the 
RT with scaling dimension A m = 2 — (2 + e/4)m 2 . There- 
fore only the lowest harmonic cos z is taken into account. 
We choose to keep 7 fixed by rewriting the renormaliza- 
tion of 7 as an additional renormalization of T, R(z) and 
v. The resulting RG flow equations, valid to lowest order 
in e, read 



IT 

dR 

H 



dv 



2-D- 2C- -v 2 A 
4 

S 2 



+ -[R"(z)] 2 -R"(0)R"(z), 



2 -2 A 



-Aw 3 . 



T = -6T, (5) 
R{z) + (zR'(z) 

(6) 
(7) 



Here I = ln(l/q) and A = -i?"(0)e 2C '. We have set 
the ultraviolet cutoff to 1 and made the substitution 
K D R(z) -> R(z) with K D X = 2 D - 1 7T D / 2 T(D/2). The 
factor exp(2(7) in the definition of A stems from the 
rescaling of the periodicity of the lattice pinning poten- 
tial. Here exp((7) has to be read as the abbreviation for 
exp(J' C,dl) with / dependent £. 

First, we remark that for D > 2 thermal fluctua- 
tions turn out to be irrelevant such that T* = at all 
fixed points of interest. For v = three locally stable 
fixed points of Eq. ([3]) are known, each with their own 
basin of attraction. These are the RF, the RB and the 
charge density wave fixed point ||(],|2^]. By inspection 
of Eq. (JsJ) it is clear that a finite v does not change 
the functional structure of the flow of R(z) since v 2 A 
approaches a finite value at the fixed point as we will 
see below. Therefore the functional form of the fixed 
point R*{z) does not depend on v and is given by a di- 
mensionless function r(u) with a characteristic scale of 
order unity. It is determined by the fixed point condition 
ar{u)+0ur'(u)+r"(u)(r"(u)/2+l) = Owithr"(0) = -I, 
where a and /3 are numerical coefficients which differ for 
RF and RB disorder, respectively. They can be deter- 
mined from the conditions that r(u) ~ |m| at large u for 
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RF's and that r(u) decays exponentially for RB's. For 
RF's we obtain a//3 = — 1, which is sufficient to deter- 
mine the exponents. For RB disorder r(0) > and we 
can choose r(0) = 1 leading to a — 1/2 by evaluating the 
above fixed point condition at u = 0. [3 — 0.6244 is deter- 
mined numerically from the condition that the function 
r(u) vanishes exponentially for u — > oo and r(0) = 1. 

To obtain the critical behavior near the RT, we as- 
sume that R(z) has already approached its functional 
fixed point form, i.e., we make the Ansatz 



R{z) = me{l)e- A «r{z/£,{iy- Ql ) 



(8) 



which is justified on sufficiently large length scales. In- 
serting this Ansatz into Eq. and using the fixed point 
condition for r(u) we obtain an equation of the form 
Ar(u) + Bur'{u) = where A and B depend on A, £ 
and its derivatives with respect to I. This equation is 
fulfilled for all u if and only if A = B = 0, which is 
equivalent to the following flow equations for A and £, 



dA 
~~dT = 



1 



2/3 , x 2 
— + 2 V 



dl 



,A 

T 



(9) 
(10) 



The rough phase of the interface is described by a sta- 
ble fixed point with v* = 0, A*/£* 2 = e/(a + 4/3). In 
this phase, RB and RF systems are characterized by a 
roughness exponent which reads in our parameter ap- 
proach £ = e/(4 + a/ (3). With the above values for a 
and P the known results ( RF = e/3 and ( RB = 0.2083e 
p0[ are exactly reproduced. For D = 4 we find that the 
interface roughness grows only sub-logarithmically with 
the system size L, (z 2 ) ~ ln CT (L/a) with a = 2/3 for 
RF's and a = 0.4166 for RB's. 

For finite v there is a new fixed point with A* = 4+e/2, 
v* = a/e/2 and £* = oo. Notice that the perturbative 
RG is still justified at this fixed point since v is small. 
Linearizing around this fixed point we obtain in lowest 
order in e the eigenvalues A± = ±2y/e. Since the fixed 
point has an unstable direction, it has to be associated 
with the RT. In the flat phase the size of typical ex- 
cursions of the interface from the preferred minimum of 
the lattice potential defines the longitudinal correlation 
length £ii ~ \v — w c | _ly " with v\\ = 1/2^/e. Interestingly 
V\\ does not dependent on a, and is therefore univer- 
sal for RF and RB disorder. The same expression for i/n 
has been obtained previously for the thermal roughening 
transition but with e = 2 — D Is) . 

The interface roughness K(x) — ([z(x) — z(0)] 2 ) can 
be obtained in terms of the effective and non-rescaled 
parameters A = Ae~ el , v = ve~ 21 using I = ln(l/q). 
Sufficiently close to the transition where £m > |x| > Lr, 
the interface roughness increases logarithmically, K(x) ~ 
ln(|x|/L fl ) on both sides of the RT, see Fig. 1(d). Be- 
yond £|| the roughness crosses over to the power law 
K(x) ~ (|x|/£||) 2 '' in the rough phase. On the flat side 
of the transition K(x) saturates on scales larger than £y 
at a finite value ~ ln(£||/L^). 



The RG calculation also confirms the result v c (Aq) = 
A(e)A^ e already mentioned above. To obtain the coef- 
ficient A(e) we integrate the flow equations (Q), (||) and 
( |To| ) by neglecting the non-linear terms in v. The ef- 
fect of the non-linear terms is subsequently included by 
matching the resulting solution v(l) and v* at the length 
scale L R , which yields A(e) = f ) 1 /C-2A( 4e )-i/C 
in good agreement with numerical solutions of the RG 
flow. Since £ < e/2, the function A(e) vanishes rapidly 
for e — > 0. Notice that in D — 4 an arbitrarily small v > 
leads to a flat phase due to the sub-logarithmic rough- 
ness for v = 0. The schematic RG flow for 2 < D < 4 is 
shown together with that expected for D < 2 and D > 4 
in Fig. 1. 

The fluctuations in the free energy grows with length 
as L e where 8 is the violation of hyperscaling exponent. 
Whereas one has exactly 9 = D — 2 + 2£ in the rough 
phase, the relevance of the lattice potential destroys the 
tilt symmetry of ([[]) leading to 9 C — 2 — e/2 at the tran- 
sition. 



D<2 





FIG. 1. Schematic RG flow as a function of the disorder 
variable A and the lattice strength v for (a) D < 2, (b) 
2 < D < 4 and (c) D > 4. Note that the fixed point as- 
sociated with the rough phase is shifted here to A = oo. (d) 
shows the length scale L dependent roughness near the RT. 

As discussed in the introduction, the RT disappears 
in the physically interesting case D — > 2+. However, the 
RT is expected to be seen even for D = 2 in systems with 
dipolar interaction. Indeed, for magnetic or ferroelectric 
domain walls with the magnetization direction p parallel 
to the wall the elastic energy has the form p3|.p4| 




(11) 



where g measures the relative strength of the dipolar in- 
teraction. Our results can also be applied to interfaces 
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with this elastic interaction. Repeating the RG analysis 
for this case, we see that the calculated exponents remain 
valid to order e if we replace e = 4— D by e = 3(3 — D)/2. 
Thus the upper critical dimension is shifted to D = 3, 
and there is a RT for 2D dipolar interfaces which is de- 
scribed by our results with e = 3/2. Note that 2D dipolar 
interfaces do not show a thermal RT [^3| . 

The threshold condition u c (Ao)Aq 2 ^ 6 = A(e) and 
hence the RT can be reached by changing the tempera- 
ture T since vq and Ajy depend in general in a different 
way on T c — T. Here T c denotes the condensation tem- 
perature of the system, e.g., the Curie temperature for 
magnets. The precise T-dependence is model dependent 
but close to T c where £o/ a 3> 1 we expect a very weak 
influence of the lattice potential, e.g., for Peierls barriers 
vo ~ e~ c ^ . Thus, for not too large disorder, we can 
expect to see this RT by increasing T. 

Under the influence of a small driving force density 
/, the interface motion is dominated by jumps between 
neighboring metastable states in the rough phase and 
between adjacent minima of the lattice potential in the 
flat phase. In both cases the creep velocity u is expo- 
nentially small @, u(f) ~ exp[-£ c /T(/ c //)' 1 ]. In the 
rough phase, » = (2( + D - 2)/(2 - C), E c ~ ^L^ 2 , 
and f c ~ 7^o-^_r 2 whereas in the fiat phase ^ = D — 1, 
E c ~ yv( 2 ~ D '' 2 and f c ~ jv for an interface with short 
range interactions. For dipolar systems, we obtain for a 
two dimensional interface in the rough phase fi RF = 1, 
l_i RB = 0.6666 and for a flat interface [i = 1 independent 
of the type of disorder. Therefore, in the RB case the 
phase of the interface can be determined by measuring 
/i. A very recent experiment on driven domain walls pH] 
shows that the exponent fi can be measured very accu- 
rately. 

To summarize, the interplay between a random po- 
tential and lattice pinning for an interface leads for 
2 < D < 4 to a continuous disorder-driven roughen- 
ing transition between a rough and a flat interface, in 
contrast to earlier results of Bouchaud and Georges for 
RB systems [JL7| which are based on a variational calcu- 
lation. However, there is a qualitative similarity in the 
phase diagram if we identify the GF phase of |l7| with 
the logarithmically rough crossover region of our results. 
Both separate the GR from the simple flat phase [f26| . 
Moreover, we have obtained the exponent vn which has 
to lowest order in e the same form as that of the thermal 
RT if e = 2 — D. For both random field and random bond 
systems, the interface shows a superuniversal logarithmic 
roughness at the transition. At the upper critical dimen- 
sion D = 4 an arbitrary small lattice pinning produces a 
flat interface due to sub-logarithmic roughness in the ab- 
sence of lattice effects. In contrast to the 2D thermal RT 
a diverging correlation length £y appears on both sides of 
the transition. On the rough side £m sets the length scale 



for the crossover from logarithmic to algebraic roughness. 
We have demonstrated that the disorder-driven RT will 
appear for two dimensional interfaces in dipolar systems. 
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Israeli Foundation (GIF) and the Volkswagen Stiftung. 
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